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Abstract. 

I review results obtained from studies of the high-redshift Lya forest, 
and present new results from HST/STIS spectra of low-redshift quasars 
in comparison with cosmological hydrodynamic simulations. The evo- 
lution of the Lya forest from z ~ 3 — > is well-described by current 
structure formation models, in which Lya forest absorbing gas at all red- 
shifts traces moderate-overdensity large-scale structures. I describe some 
of the insights provided by hydrodynamic simulations into the observed 
statistical trends of Lya absorbers and the physical state of the absorbing 
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. 1. Introduction 
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Quasar absorption lines have been used as probes of the high-redshift universe 
for some time (e.g. Sargent et al. 1980), since the Lya (1216A) transition is 
redshifted into the optical at z £ 2. HIRES (Vogt et al. 1994) on the Keck 10m 
telescope has been at the forefront in providing high S/N, high resolution quasar 
spectra that has enabled dramatic advances in our understanding of the high- 
redshift intergalactic medium (Rauch 1998). Conversely, studies of low-redshift 
Lya absorbers require ultraviolet spectroscopy, which has only recently become 
routinely possible with the Hubble Space Telescope. The HST Quasar Aborption 
Line Key Project (Bahcall et al. 1993) observed over 80 quasars using the Faint 
Object Spectrograph (FOS), and found a surprisingly large number of absorbers 
as compared to an extrapolation from high-redshift (Bahcall et al. 1996; Wey- 
mann et al. 1998). Recently the deployment of the Space Telescope Imaging 
Spectrograph (STIS) has provided a dramatic increase in resolving power (res- 
olution 7 kms" 1 as opposed to FOS's 230 kms -1 ), thereby fully resolving all 
Lya absorbers seen in UV quasar spectra. While only the brightest quasars 
and AGN may be easily observed with STIS, this has still resulted in significant 
gains in our understanding of low-redshift Lya absorbers and their relation to 
high-redshift absorbers. In these proceedings I report some recent results from 
HST/STIS. 

In conjunction with these observations, hydrodynamic simulations of struc- 
ture formation in currently favored cosmological scenarios have elucidated a 
new paradigm for the nature of Lya forest absorbers, particularly at high red- 
shift. These simulations indicate that Lya absorbers arise in highly photoionized 
diffuse intergalactic gas tracing the dark matter in non-equilibrium large-scale 
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structures. Their temperature is set by a balance between photoioinization heat- 
ing from a metagalactic UV background (presumably from quasars; Haardt Sz 
Madau 1996) and adiabatic cooling due to Hubble expansion, resulting in a 
tight relation T oc p 0,6 . These simulations are able to reproduce various statis- 
tical properties of high-redshift Lya absorbers in detail, spanning the range of 
damped Lya systems to the weakest observed detected with HIRES, all within a 
model that generically arises in CDM cosmologies (Hernquist et al. 1996; Dave 
et al. 1997). Figure 1 shows a comparison of a portion of the HIRES spectrum 
of Q1422+231 (Songaila & Cowie 1996) with a simulated spectrum covering the 
same redshift interval, having resolution and noise properties that emulate the 
data. Can you tell which is which? (answer at the end). 

Due to the relatively simple relationship between H I and the underlying 
dark matter, a formula known as the "Fluctuating Gunn-Peterson Approxima- 
tion" (FGPA; see Croft et al. 1998 for full formula) provides a remarkably 
accurate description of Lya absorbing gas: 

thi oc p^Tm' 1 , (1) 

where thi is the H I optical depth, p is the density of dark matter (or baryons), 
and Thi is the metagalactic H I photoionization rate incident on the absorber. 

2. Results From High-Redshift Lya Forest Studies 

The FGPA has been combined with HIRES quasar spectra to provide stringent 
constraints on many aspects of the physics of high-redshift intergalactic gas. For 
instance, a measurement of the mean optical depth together with an independent 
estimate of Thi result in a robust estimate of ss 0.02/i~ 2 (Rauch et al. 1997), 
in good agreement with D/H measurements (Buries & Tytler 1998) and recent 
CMB analyses (e.g. Pryke et al. 2001). 

The FGPA also implies that a given Lya optical depth (or flux, in the 
optically-thin regime) corresponds to a particular density of gas at a particular 
temperature. Hence ionization corrections can be accurately obtained to con- 
strain the metallicity of the diffuse high-redshift intergalactic medium (IGM) 
from metal-line observations of C IV (Songaila & Cowie 1996; Rauch, Haehnelt 
& Steinmetz 1997; Dave et al. 1998). These studies found [C/H] ss -2.5 for a 
Haardt & Madau (1996) ionizing background shape, and significantly lower for 
a spectrum with fewer high-energy photons (as would be expected if Helium had 
not yet reionized; see Heap et al. 2000). Furthermore, recent VLT/UVES ob- 
servations of O VI absorption suggests the presence of some metals even in voids 
(Schaye et al. 2000), presenting a significant theoretical challenge to models of 
metal ejection and transport in the IGM (e.g. Aguirre et al. 2001). 

The most ambitious use of the FGPA, and the one for which it was origi- 
nally developed, is a reconstruction of the matter power spectrum at z ~ 3 from 
Keck/HIRES quasar spectra. Since the optical depth in the Lya forest traces the 
matter density, fluctuations in the optical depth reflect a 1-D probe of the am- 
plitude and shape of the matter power spectrum, at scales of ~ 1 — > 10/i _1 Mpc 
(comoving). For a flat universe, this measurement yields a constraint on the 
matter density of Q m ~ 0.5 + 0.29(r — 0.15), where T is the power spectrum 
shape parameter (Croft et al. 2000). This value is slightly higher that in the 
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"concordance model" (Q m ss 1/3; Bahcall et al. 1999), but the latest CMB 
measurements combined with the Hq Key Project value (Mould et al. 2000) 
suggest a higher Q m as well (Pryke et al. 2001). Thus the Lya forest provides 
an independent avenue for doing "precision cosmology", alongside CMB, Type 
la supernovae, etc. 




Figure 1. 100A section from the spectrum of Q1422+231, and an 
artificial spectrum drawn from a hydrodynamic simulation of a ACDM 
universe, having resolution and noise characteristics constructed to be 
similar to Q1422. The fluctuating absorption pattern seen in the sim- 
ulated spectrum (which one is it?) arises naturally from large-scale 
structure generated by hierachical collapse, and is statistically indis- 
tinguishable from that seen in Q1422. 
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3. Evolution of the IGM to Low Redshift 

Recently, hydrodynamic simulations have been evolved to redshift zero to in- 
vestigate the nature of low-redshift Lya absorbers and the evolution of the Lya 
forest. They indicate that while many of the baryons have moved into galax- 
ies, clusters, and a diffuse shocked "warm-hot intergalactic medium" (Cen & 
Ostriker 1999; Dave et al. 2001), roughly one-third of all baryons continue to 
reside in photoionized gas tracing large-scale structure. 

These simulations are able to explain a wide range of Key Project obser- 
vations quite naturally within the context of hierarchical structure formation 
scenarios. For instance, the Key Project data showed that high-redshift (z ^ 2) 
Lya absorbers are evolving away quite rapidly, whereas for z <; 1.5 the absorber 
number density evolution slows abruptly and dramatically (Bahcall et al. 1996). 
In the simulations of Dave et al. (1999) and Theuns et al. (1998), this occurs 
due to the dimunition of the quasar population providing the metagalactic pho- 
toionizing flux. As a result, the rapid evolution of absorbers due to Hubble 
expansion in the high-z universe is countered by the increase in neutral fraction 
of Lya absorbing gas at z <J 2. These simulations disfavor a scenario in which a 
different population of absorbers dominates at low redshift as compared to high 
redshift, and instead suggest that the nature of absorbing gas is quite similar. 




Figure 2. Top panels: Relation between column density of absorber 
and maximum physical density (in units of mean cosmic density) at 
the absorber peak optical depth, at z = 3 and z = 0. Best-fit relation 
is shown as the solid line, and is reproduced in the other panel as the 
dotted line. Bottom panel: For three selected densities, the evolution 
of the corresponding absorber column density as a function of redshift. 
The z ~ forest is qualitatively quite similar to the z ~ 3 forest, except 
column densities are shifted down by ~ 1.5 dex. 
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Weymann et al. (1998) found that stronger absorbers disappear more 
rapidly than weaker absorbers, and in fact the weakest absorbers seen in the 
FOS sample (rest equivalent width W r ph 0.1 A) show an increase in number 
between z ~ 1.5 — » 0. Dave et al. suggested that this arose due to differen- 
tial Hubble expansion, in particular that low-density regions expand faster than 
high-density ones and therefore provide an increasing fractional cross-section of 
absorption to lower redshifts. 

Motivated by the good agreement between simulations and FOS observa- 
tions, Dave et al. investigated the physical state of the gas giving rise to Lya 
absorption at low redshift. They found that the gas follows similar physical 
relationships as deduced at high-redshift, and is thus of the same basic charac- 
ter, but that a given column density absorber corresponds to different physical 
densities at different redshifts. 

This effect is shown in Figure 2. The top left panel shows the density- 
column density relation at z = 3 from a ACDM simulation. The tight relation 
reflects the accuracy of the FGPA; much of the scatter about the best-fit relation 
(solid line) comes from errors in Voigt profile fitting (noise with S/N=30 was 
added to the artificial spectra to roughly match observations). At low redshift, 
the relationship persists, though there is more scatter due to the presence of 
absorbers arising in shock-heated intergalactic gas (Dave et al. 1999; Cen & 
Ostriker 1999). The best-fit slope of the p — Nm relation is similar, but the 
amplitude has shifted considerably from z = 3 — > 0. This shift arises because 
of the interplay between Hubble expansion and the evolution of the ionizing 
background. 

The bottom panel of Figure 2 focusses on the evolution of this trend with 
redshift. Here, we show the typical column density of an absorber arising in 
gas with three different densities, as a function of redshift. For instance, gas at 
the mean density (which typically remains around that density over a Hubble 
time) produces a strong absorber of Nm 10 14 cm -2 atr z = 3, but by z = 
it produces only a very weak absorber of Nm ss 10 12 ' 5 cm~ 2 . In general, a given 
density corresponds to a column density ~ x30 lower at z = as compared to 
z = 3. Thus studying physically and dynamically equivalent absorbers at high 
and low redshift requires comparing across different column densities. 

Using this trend, the knowledge and intuition gained from studies of the 
Lya forest at high redshift can be translated to expectations at low redshift. 
For instance, it is known that most absorbers with A%i ;> 10 145 cm -2 at z ~ 3 
are enriched (Songaila & Cowie 1996); thus one expects that absorbers with 
Nm ;> 10 13 cm~ 2 should be enriched to at least a similar level (though detecting 
these metal lines will be challenging). 

4. Low-Redshift Lya Forest Observed with HST/STIS 

A new window on the low-redshift Lya forest has opened with the deployment 
of STIS aboard HST. STIS's E140M Echelle grating provides ~ 7 kms" 1 reso- 
lution, a level of detail comparable to HIRES albeit with somewhat more noise. 
STIS can thus probe weak forest absorbers that are predicted to be physically 
similar to those at high-redshift, and allows a detailed comparison to be made 
with hydrodynamic simulations. 
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Dave k, Tripp (2001) compared the quasar spectra of PG0953+415 and 
H1821+643 to carefully-constructed artificial spectra drawn from a ACDM sim- 
ulation. The statistical properties of absorbers in the simulations were in re- 
markable agreement with observations, as shown in Figure 3. The column den- 
sity distributions agree down to the smallest observable absorbers, in both slope 
and amplitude. Interestingly, the slope of the column density distribution (from 

dzdNm ^ ^ui) ls measured to be (3 « 2.0 ± 0.2, which is considerably steeper 
than at high redshift ((3 « 1.5; Kim et al. 2000). This indicates, in agreement 
with Weymann et al., that stronger absorbers have evolved away faster than 
weaker ones. Our result is also in broad agreement (~ la higher) with Penton 
et al. (2000), who found (3 » 1.8 from GHRS data. 
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Figure 3. Column density (left panel) and 6-parameter (right panel) 
distributions from HST/STIS observations of PG0953 and H1821 (solid 
lines), compared with artificial spectra (dashed lines). Dotted line in 
left panel shows the result of increasing S /N by x 5 in the artificial spec- 
tra, indicating absorber sample is complete down to Nm <J 10 13 cm -2 . 
Dotted line in right panel shows the simulated line widths purely from 
thermal broadening, indicating that thermal broadening provides a 
substantial component of line widths at low redshift. 

The distribution of line widths (6-parameters) is also in excellent agreement 
with observations. The median 6-parameter is around 21 kms -1 , but this value is 
sensitive to the column density detection threshold, because the column density 
and 6-parameter are correlated. This can be understood from the fact that higher 
column density systems arise in higher density gas that has higher temperature. 
Since thermal broadening contributes significantly to line widths (compare the 
dotted and dashed lines in the right panel of Figure 3), higher column density 
absorbers will tend to have larger widths. Thus Shull et al. (2000), who find 
a median 6-parameter of « 30 kms -1 from FUSE and GHRS data, are not in 
disagreement with our result because their sample has a higher column density 
detection limit. 

The large thermal component of linewidths facilitates a measurement of 
the typical temperature of Lya absorbing gas at low redshift. By comparing the 
observed median 6-parameter to simulations, Dave & Tripp found that purely 
photoionized Lya absorbers (those that have not been shock-heated) at the mean 
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density have a temperature of T ~ 5000K, in agreement with Ricotti, Gnedin & 
Shull (2000). 

The agreement in amplitude of the column density distribution can be trans- 
lated into a constraint on the H I photoionization rate Thi, since varying Thi 
would shift the simulated column density distribution proportionally in the hor- 
izontal direction (assuming optically thin lines, which is valid for most STIS 
absorbers). Dave & Tripp obtain Thi ~ jq-13.3±o.7 g -i^ w i^jj errors domi- 
nated by systematic uncertainties in modeling, providing among the most precise 
constraints on the metagalactic Thi to date. 

5. Conclusion 

The Lya forest has proved to be a powerful tool for constraining a wide variety 
of physical and cosmological parameters, when hydrodynamic simulations and 
the FGPA are used in conjunction with high-resolution quasar spectra. Much 
progress has already been made at high redshift, where the metallicity, ionization 
state, and mass fluctuation spectrum have all been derived from Lya forest data. 

Hydrodynamic simulations of structure formation also provide an excellent 
description of low redshift Lya forest absorbers, particularly for the weak ab- 
sorbers that reside in moderate-overdensity regions analogous to high-redshift 
absorers. HST/STIS has opened up a new era in the study of the low-redshift 
IGM, which promises to progress considerably further with the deployment of 
the Cosmic Origins Spectrogaph aboard HST (in 2003). Initial results indicate 
that the low-redshift Lya forest can yield important constraints on physical pa- 
rameters in the IGM, as well as provide a more complete picture of the evolution 
of baryons in the Universe. 
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